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The missing compact star of SN1987A: a solid quark star? 
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\ 1 INTRODUCTION 



ABSTRACT 

To investigate the missing compact star of Supernova 1987 A, we analyzed both the cooling 
and the heating processes of a possible compact star based on the upper limit of observational 
X-ray luminosity. From the cooling process we found that a solid quark-cluster star, which has 
a stiffer equation of state than that of conventional liquid quark star, has a heat capacity much 
smaller than a neutron star. It can cool down quickly, which can naturally explain the non- 
detection of a point source (neutron star or quark star) in X-ray band. On the other hand, we 
consider the heating process from magnetospheric activity and possible accretion, and obtain 
some constraints to the parameters of a possible pulsar. We conclude that a solid quark-cluster 
star can be fine with the observational limit in a large and acceptable parameter space. A 
pulsar with a short period and a strong magnetic field (or with a long period and a weak field) 
would has luminosity higher than the luminosity limit if the optical depth is not large enough 
to hide the compact star. The constraints of the pulsar parameters can be tested if the central 
compact object in 1987 A is discovered by advanced facilities in the future. 
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Pulsars are thought to be neutron stars (NSs) since the discovery 
of the first pulsar dHewishetal1 [l968). Nevertheless, th e equa- 
tion o f state (EoS) of neutron stars is not clear till now. IWittenl 
J 19841) pointed out that the true ground state of hadrons may be 
strange matter, which contains roughly equal numbers of up, down 
and strange quarks. It is then realized that a quark star or strange 
star (SS) might be the ground state of NSs dAlcock et ail 1 19861 : 
lHaensel et aflll986h . NSs may convert to SSs. 

Are pulsars NSs or SSs? Essentially, this is a problem of non- 
perturbative quantum chromodynamics (QCD), which describes 
the strong interaction processes in low energy scales, and is very 
hard to solve in mathematics. Fortunately, the astronomical obser- 
vations could discriminate NSs and SSs, and help to understand the 
non-perturbative QCD issue. 

There are many methods to constrain the EoSs of pulsars and 
check the differences between NSs and SSs. For example, SSs are 
self-constrainted, while NSs are gravity constrainted, so that they 
have different mass-radius relations. If we can measure the mass 
and radius of a pulsar, we can discriminate whether it is a NS or a 
SS. This has not succeeded up to now, mainly because of uncertain- 
ties of the radii. Nevertheless. iLi et alj dl999h suggested that SAX 
11808.4-3658 may be a quark star based on this relation and avail- 
able data. Another constraint is the minimum mass of a pulsar or 
sub-millisecond pulsar. Since SSs are self-bounded, their mass can 
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be much smaller than the minimum mass of NS, which could be 
about O.IMq and can s pin with a per iod shorter than 1 millisec- 
ond while a NS can not jDu et al.1l2009h - If the mass of a pulsar is 
smaller than 0.1 Mq or the period is shorter than 1 millisecond, it 
could not be a NS. 

Different EoS of a pulsars gives a different maximum 
mass. The newly con firmed mass of PSR 11614-2230 of 2Mq 
jPemorest et alj |2010) rules out al most all currently propos ed hy- 
peron or boson condensate EoSs ( Lattimer & Prakash 2007) and 
traditional soft EoS of SS dChan et al.|[2009h . lLai & Xul j2009t) ar- 
gued that the solid quark-cluster stars (SQS), a special kind of SS, 
could have the maximum mass to be larger than 2 Mq because of 
the very stiff EoS. 

The bolometric radiation of a young pulsar can be used to dis- 
tinguish NSs and SSs, because they are different in thermal capacity 
and surface radiation. We note that a SQS has very small thermal 
capacity because its temp erature is much lower than the Debye tem- 
perature dYu & Xii201ll) . it can cool down very fast and has a very 
lower bolometric luminosity than a NS. The missing compact ob- 
ject in SN1987A dMcCravl l2007) is studied in this paper with this 
luminosity constraint. 

There is no doubt that SN1987A provides an unprecedented 
opportunity in astronomy and astrophysics studies. Yet the pos- 
sible compact star produced in this explosion is still a mystery. 
On 23 February 1987, both the Kamiokande II detector and the 
Irvine-Michigan-Brookhaven detector observed a neutrino burst 
dHirata et al]|l987l : iBionta et"al] [l987) from the Large Magellanic 
Cloud, just before the optical shine. The energy spectrum and the 
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flux density of the neutrino burst are consistent with the current the- 
ory of core-collapse supernova with an energy release of ^3x 10 53 
ergs, from which a NS is expected to form therein. Astronomers 
were excited by the possibility of watching the NS birth at the very 
beginning. From then on, the powerful detectors or telescopes on 
the ground and in the space searched it from radio to gamma ray 
bands, one after another. Unfo rtunately, no pulses were detected or 
no point source is found there (McCrav 2007; Manc hestej|2007h . 

Parkes 64-m radio telescope was used to search pulses from 
SN1987A in the first few years at frequencies betw een 400 and 
5000 MHz, down to the limit of ~0.2 mjy at 1.5 GHz (Man chester! 
The strong efforts were made with Parkes in 2006 in sev- 



eral bands, but no pulsar ca ndidates with a S/N ratio greater than 
9.0 were found (Manc hesteJl2007h . Some opti cal observations an- 
nounced the detection of a pulsar in SN1987A dKristian et al . 1989; 
iMurdinl 1 199dl; iMiddleditchl |2000h . but all of them have not been 
confirmed. iPercival et alj d 1995T) made an optical search using the 
Hubble Space Telescope. No significant pulsations were found in 
the period range of 0.2 ms to 10 s with an upper limit for the 
pulsed emission equi valent to a V magnitude of ~24. A similar 
search was made by Manch ester & Petersonl J 19961) using the 3.9-m 
Anglo- Aust ralian Telescope w i th sim ilar parameters, and they got 
null results. Shtvkov skiv et alj |2005) obtained a luminosity limit 
in the 2-10 keV band o f 5 x 10 34 erg s" 1 using the XMM-Newton. 
IPark et alj l2002ll2004f) obtained upper limits of 5.5 x 10 33 erg s" 1 
and 1.5 x 10 34 erg s _1 in the same X-ray band using Chandra. 

To explain the non-detection of the predicted NS, several pos- 
sibilities were discussed. First, some of the ejected material of SN 
may fall back to the NS surface shortly after the supernova explo- 
sion. We can not see the NS since it has been converted to a BH 
due to the fall back material. Bu t this possibility sho uld be small. 
The discovery of a 2 Mq pulsar dDemorest et alj20 10) suggests the 
maximum mass of pulsar could be much larger than 1 .4 Mq , and 
therefore a normal NS needs to accrete a lot of material for the con- 
version to a BH. Anoth er possibility is that the NS may be located 
in the cold dust cloud dMcCravl 12007) at the center of SN1987A, 
which may be opaque in some bands. But if the NS is not inside 
the dust or the dust is optical thin in X-ray band, it could be very 
intriguing because the observed upper limit at around 1 keV is al- 
ready lo wer than the expected luminosity of a cooling NS w ith no 
heating dPark et alj|2007l ; iMcCravl 120071) . Ichan et all |2009|) sug- 
gested that the compact star at the center of SN1987A may be not 
a NS but a SS with a softer EoS than that of neutron star matter, 
which could have a X-ray luminosity less than 10 34 erg s -1 at the 
age of 20 years. Note that the S S of this EoS is alm ost ruled out 
by the newly found 2 Mq NS dPemorest et al.l l2010). because its 
maximum mass is smaller than 2 Mq due to its soft EoS. 

We emphasize that a model for quark star with a stiff EoS, e.g. 
the SQS jLai&Xdl201lh. is still survived the oretically. The SQ S 
is suggested by Ixd J2003h and impro ved by IZhou et al] d2004) . 
lLai & Xj d2009T) and lLai & Xul ( 1201 lh . The key point is that the 
interaction energy of quarks in a compact star could be higher than 
the Fermi energy when t he density is lower than a few tens of nu- 
clear densities ( Xu 201 1), so that quarks may be clustered and the 
star could be a solid star, i.e. the solid quark star (or quark-cluster 
star). The SQS model can naturally explain most of the observa- 
tional features of pulsars dXull201 If) . The surface of SQS is self- 
bounded which provides a larger binding energy than the gravity- 
bounded NS. That is more hel pful to the gener ation of sparks for 
radio emission dXu et al. 1999; Oiao et al. 2004). The X-ray emis- 
sion and pulsations in magnetosphere are similar to those of a NS. 
The star quake of SQS can induce glitches and energetic bursts 



which may be shown as soft gamma-ray repeaters dXu et al]2006T) . 
The SQS co uld has a larger maximum mass than other model of 
quark stars. lLai & Xul d2009h gives a ma ximum mass > 2Mq , 
which stands the test of the 2 Mq pulsar dPemorest et alJr2010h . 
The SQS has a smaller radius than any NSs when the mass is 
small. The newly discover y of the radius and mass of the rapid 
burster (MXB 1730-335) dSala et alj[2012h could naturally fit the 
SQS model. Ad ditionally, the SQS has a very low heat capacity 
{Si & Xvj[2oTlh which helps to cool down very fast. This may ex- 
plain the non-detection of a point source in SN1987A. 

In this work we will first analyze the cooling process of a SQS 
in comparison with NS and traditional liquid SS, and then further 
study the constraints on the parameters via heating processes of 
pulsars. Our conclusions and discussions are presented in the last 
section. 



2 COOLING OF THE POSSIBLE COMPACT OBJECT 

A stellar BH has no classical radiation, and its Hawking radiation is 
negligible in astronomy, thus its cooling luminosity is nearly zero. 
For NSs and SSs, the cooling processes are determined by their 
heat capacities and surface radiations. 

The heat capacities of NSs, conventiona l SSs and SQSs are 
dMaxwelllll97^ ; lNg et al.ll2003l ; lYu & Xvj|201ll) . 



Cns — C N s + Cns, 

C SS = C s q s + Cfs 7 + C C 3S , 

Csqs = Csqq + Csqsj 



(1) 
(2) 
(3) 



where the superscripts n, e, q, g-7, and 1 denote the contributions of 
neutrons, electrons, quarks, quark-gluon plasma and lattice struc- 
ture, respectively. In eqs. (1) and (2), C NS an d Cg S are larger than 
C NS an d Css when temperature is higher than critical temperature 
T c (~10 9 K). When T < T c the superfluid state appears, C NS 
and Cg S will exponential decay and vanish quickly. In eqs. (2) and 
(3), both C|g~ 7 and C S q S are in proportion to T 3 , while Cgg and 
Cgqg are in proportion to T. Thus the heat capacity of compact 
stars is dominated by electrons when temperature is not too high, 
with T c ~10 9 K for NS and SS, and ~10 10 K for SQS dYu & Xul 
l]). The temperature will otherwise quickly cool down below 
10 9 K within dozens of seconds. Therefore, the heat capacity of 
electrons is overwhelming the cooling process in almost all obser- 
vati onal time. The heat capacity of electrons in N S, SS and SQS 
are dMaxwelllll979l ; lNg et al.ll2003l;lYu & Xull201lh . 



C£ s = 1.9 x 10 37 M lP ^ 3 T 9 erg K"\ 

cgg = 1.7 x W 20 (Yp/p ) 2/3 T 9 erg (cm 3 K)" 1 , 



C 



SQS 



k B T s 



(4) 
(5) 

(6) 



where Mi = M/Mq, p u = p/10 14 g cm -3 , T 9 = T/10 9 K, 
Y is the electron fraction, po is nuclear matter density, N e is the 
electron number in a star, fcs is Boltzmann's constant, T s is the 
value in the star's local reference frame, and Ef is the Fermi energy 
of the degenerate electron gas. In the extremely relativistic case, 
Ep = ( 3n g^ ) 1,/3 c, where n e is the number density of electrons, 
h is the Planck constant, and c is the speed of light. From eq. (6) 
we have 



Csqs * 3.5 x 10 37 (YM^ 3 R ti T 9 erg KT 1 , 



(7) 
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Figure 1. The cooling curves of SQSs with bremsstrahlung (BR) 
ICaron & Zhitnitskvl 120091) and black body radiation (BB). Neutrino radi- 
ation and color superconductivity related photon emission are not consid- 
ered. The observational upper limit is indicated as a horizontal dotted line. 
It shows that the cooling luminosity of a SQS could be smaller than 10 34 
ergs s — 1 about 20 years after its birth even it cools down by bremsstrahlung 
emission. Here we take the stellar mass M = Mi Mq and the number ratio 
of electron to baryon as Y . 



where Re is the star radius in units of 10 cm. The electron number 
in quark star is much smaller than that in neutron star, Y ~ 10~ . 
Thus the heat energy conserved in a quark star is far less than in a 
neutron star. 

Surface radiation is important in the cooling process. On the 
NS surface, the radiation can be simply treated as approximate 
black body radiation because of atomic atmosphere. On SS or 
SQS surface there is no atomic atmosphere, therefore the radia- 
tion mai nly depends on the interaction between electric layer and 
photons. IChanetal.1 b009l) suggested that the surface radiation 
of SS is via bremsstrahlung, and predicted a luminosity smaller 
than 10 34 erg s _1 when it is older than 20 years. In the ir pa- 
per the bremsstrahlun g calculations o f Jaikumar et al.l J2004l) were 
used. We noticed that ICaron & Zhitnitskvl 1(2009!) got a much low 
flux of the S S bremsstrahlung emission . In our work we adopt 
the formula of ICaron & Zhitnitskvl d2009h to calculate the cooling 
of SQS, which makes our results more reliable. Figure 1 shows 
the bremsstrahlung cooling curves of S QSs toge t her w ith those 
of black body radiation of the model of lLai et all feOllh . We ig- 
nored neutrino radiation and color superconductivity related pho- 
ton emission mechanisms. It shows that the cooling luminosity of 
a SQS could be smaller than 10 34 ergs s _1 about 20 years after its 
birth even it cools down by bremsstrahlung emission. The SQS was 
not detected previously, when it has a larger luminosity, because of 
the deep X-ray optical depth dMcCray| 2007) and poor observation 
devices. Thus we can conclude that if there is a cooling SQS in 
SN1987A it should be undetectable up to now, which is compatible 
with the observations. 



3 CONSTRAINTS ON THE PARAMETERS VIA 
HEATING PROCESSES 

The bolometric luminosity of a compact star comes not only the 
contribution of cooling process but also heating processes from 
the activity of magnetosphere and the accretion from interstellar 



medium and accretion disks. The heating luminosity is probably 
lower than the upper limit of observations. The heating mechanisms 
are independent of the EoS of pulsar, therefore we can make some 
constraints on the physical parameters of a possible pulsar regard- 
less it is a NS or a SQS. 

Usually, the activity of magnetosphere can produce both ther- 
mal and non-thermal X-ray emissions. When particles are acceler- 
ated in the magnetosphere they can emit non-thermal X-rays with 
a power law spectrum. When particles flow in and bombard on 
the pulsar surface, t hey will heat it, and therm al X-rays come out 
from the heated area. iBecker & TriimpeiU 19791) found that the non- 
thermal X-ray luminosity L x and spi n down energy loss rate E 
are roughly related by L x = 10~ 3 g. lYu & Xul feOllh found that 
the thermal X-ray luminosity has a similar relation with that E, 
LZi ~ io _3 £. Since the observations have give n a upper limit of 
therm al X-ray luminosity LZi ~ 10 34 erg s" 1 JPark et alj|2002l 
2004), we can get an upper limit of spin energy loss rate E. And 
we can further make a constraint to the parameters of spin and mag- 
netic field of the possible pulsar, because E is the function of spin 
and magnetic field, 

2 



E 



3c 3 



A (V 1 



(8) 



where n± and are the vertical fraction of magnetic moment /i 
and the spin angular frequency of the pulsar. We use a parameter 
a for the factor of thermal X-ray luminosity in terms of spin down 
energy lost rate, i.e. 



And then we get 



-1/2 ; 



1/2, 



(9) 



(10) 



where P = 2n/Q is the spin period of a pulsar. 

Note that most of pulsars have an X-ray luminosity in the 
region of 10~ 4 — 10~ 2 E in both therma l and nonthermal cases 
jBecker & Trumped 1 19791 : lYu & Xdfcoflh . We used a = 10" 4 , 
10 -3 and 10 -2 to get the upper left three lines in two panels of 
Figure 2. In the figure, the region to the left of the dash-dotted line 
should be ruled out because even the X-ray factor is as small as 
10 -4 , it's luminosity should be larger than 10 34 erg s _1 . The re- 
gion below the solid line should be compatible with observational 
upper limit, because the X-ray luminosity is smaller than 10 34 erg 
s _1 even a = 10~ 2 . 

The accretion heating can also be used to restrict pulsar pa- 
rameters. The accretion from interstellar medium is intensively de- 
pendent on the proper motion of a pulsar, and usually has a very 
small luminosity which we can ignore. Here we consider only the 
possible accretion from the fall back disk. 

To understand the evolution of the fall back disk we define 
three radius: the light speed radii rt, the co-rotation radius r co and 
the magnetosphere radius r m , as, 



r L = cP/2tt, 

r co = (GM/4tt 2 ) 1/3 P 2/3 , 
r m = fi 4/7 (2GM)- 1/7 M- 2/r , 



(11) 
(12) 
(13) 



where G, M and M are the gravitational constant, the mass of pul- 
sar and the accretion rate, respectively. In fall back disk regime a 
pulsar often undergoes three phases: pulse phase, propeller phase 
and accretion phase. The pulse phase happens in the early stage of 
a pulsar, when a young pulsar has a fast spin and strong radiation 
which pushes the interstellar medium out of the light speed radius 
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Figure 2. Constraints on parameters of the possible compact object in 
SN1987A via magnetosphere action heating (upper left three lines) and fall 
back disk accretion heating (lower right three lines). Left panel: Mi = 2 
and b = 0.9, a = 10~ 4 for red dash-dotted line, a = 10 -3 for red dashed 
line, a = 10~ 2 for red solid line, M = 10 18 erg s" 1 for blue solid line, 
M = 10 16 erg s _1 for blue dashed line and M = 10 14 erg s — 1 for blue 
dash-dotted line. Right panel: Same as the left but Mi = 1. 



r_L, thus the compact star acts as a pulsar. When the radiation be- 
comes not so strong, some of the medium may go into the light 
radius and interact with the magnetosphere. When r co < r m < tl 
or even r m is a little smaller than r co , the magnetic freezing effect 
would force the flow in medium to co-rotate with the pulsar, and the 
spin energy of a pulsar will transfer to the medium via the interac- 
tion. This is so called the propeller phase. In this phase only small 
amount of medium can diffuse and fall on the star surface, which 
induces the low luminosity X-ray emission. If the accretion mate- 
rial can across the co-rotation radius and access the so called break 
radius r or , which is a little smaller than r co (i.e. r m < 7V < r co ), 
a massive accretion would occur. This is the accretion phase, which 
usually produces a very large X-ray luminosity. 

The Hubble Spac e Telescope f ound dust clouds interior to the 
debris in SN1987A dMcCravl 120071) . which could be indication of 
the formation of a fall back disk around the centra compact star. If 
the fall back disk does exist, we would expect the pulsar is not in 
accretion phase which otherwise results in X-ray luminosity larger 
than 10 34 erg s _1 . Thus the disk should have r m > rv. We pre- 
sume 7V is proportion to r co , i.e. r or = br co , here b is a constant 
and < b < 1. Therefore to fit with the observations we need 



> br c 



(14) 



Submitting r m and r co to the above inequation, we get the relation 
of the pulsar parameters, 

Mao > 0.0746 7/4 A/ 1 5/6 A/ 1 1 6 /2 P 7/6 , (15) 

where /x 30 = ft/W 30 and M 16 = M/10 16 . 

The lower right three lines in both panels of Figure 2 show 
the magnetic moment depend on the spin period with other given 
parameters. We presumed b — 0.9 in both cases, Mi = 2 for left 
and Mi = 1 for right. We used the accretion rate from M = 10 18 g 
s _1 to M = 10 14 g s^ 1 . It should be noticed that when M < 10 14 
g s _1 the heating luminosity can not exceed 10 34 erg s _1 even in 
the accretion phase. For a disk around a pulsar the accretion rate is 
almost impossible to exceed 10 18 erg s^ 1 . 

Therefore, two kinds of heating mechanisms provide two sets 



of limits to the parameters of a possible pulsar, as shown in Figure 
2. From the combination of both heating mechanisms, the region 
between two solid lines is fine with observations, while the upper 
left and lower right regions are the forbidden areas. From the com- 
parison of left panel and right panel we can see that pulsars with 
a smaller mass has a wider region of parameters than the massive 
one. Because the heating mechanisms are EoS independent, our re- 
sults are applicable to both NS and SQS. 



4 CONCLUSIONS AND DISCUSSIONS 

We analyzed the cooling processes of compact stars and found that 
the SQSs could have a thermal luminosity lower than the obser- 
vational upper limit of SN1987A. We studied the possible heating 
processes of a young pulsar and obtained some constraints on the 
pulsar parameters. It should be reasonable to conclude that: 

(i) A solid quark-cluster star with normal parameters is compat- 
ible with the non-detection in SN1987A, because both its cooling 
and heating luminosities should be lower than the observational up- 
per limit. A low mass quark star has a wider parameter space than 
a massive one. 

(ii) If the compact star is shielded by dust, the parameter con- 
straints should be relaxed. The parameter space of a SQS can be 
wider, even the object could be a normal NS. 

(iii) However, a black hole candidate can not be ruled out by 
cooling and heating analysis, though it is very difficult to form via 
accretion as we mentioned before. 

The predicted parameter space in this paper could be tested 
by the future observations. The impact of the supernova blast wave 
with its circumstellar matter is producing a ring which is visible 
from mm-band to X - ray band of SN1987A jBouchet et alj|200rj : 
lGaensleretaIll2007l ; iMcCravl [20071 : iFrance et all 1201 d) and it is 
still brightening. This ring would make it almost impossible to de- 
tect or rule out a cooling compact star in these bands. But a heating 
pulsar may be detected in the future because the heating luminosity 
from magnetospheric activity decays very slowly, and when the fall 
back disk evolves into the accretion phase so that the X-ray lumi- 
nosity could be much larger than that in other phases and finally it 
could be detectable. It is not clear whether the low frequency radio 
and high energy 7-ray bands are affected by the ring which pre- 
serves the possibility of discovering the compact star and test the 
parameters by the future advanced facilities, e.g. the Square Kilo- 
meter Array (SKA) telescope. 
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